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ABSTRACT 

We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and 
its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the 
HI 21 cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS 
Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at 
impact parameters of 63 to 231 kpc. The spectra were obtained with the Cosmic Origins Spectrograph 
aboard the Hubble Space Telescope. We detected neutral hydrogen (Lya absorption-lines) in the CGM 
of 92% of the galaxies. We find the radial profile of the CGM as traced by the Lya equivalent width 
can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. 
We found no correlation between the orientation of sightline relative to the galaxy major axis and the 
Lya equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the 
atomic gas in the interstellar medium. We find a strong correlation (99.8% confidence) between the gas 
fraction (M(HI)/M*) and the impact-parameter-corrected Lya equivalent width. This is stronger than 
the analogous correlation between corrected Lya equivalent width and SFR/M* (97.5% confidence). 
These results imply a physical connection between the H I disk and the CGM, which is on scales an 
order-of-magnitude larger. This is consistent with the picture in which the H I disk is nourished by 
accretion of gas from the CGM. 

Subject headings: galaxies: halos — galaxies: ISM — quasars: absorption lines 


1. INTRODUCTION In the standard paradigm, galaxies grow primarily 

through the accretion of gas that flows from the Inter- 
Galactic Medium (IGM). Much of this accreted gas is ul¬ 
timately turned into new stars. In turn, the massive stars 
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release matter and energy that can affect not only the 
galaxy itself, but also the IGM. These flows in and out 
of galaxies pass through the Circum-Galactic Medium 
(GGM), a region extending out from the galaxy itself to 
roughly the virial radius of the surrounding dark matter 
halo. Th e GGM also repre s ents a significant res ervoir of 
baryons (|Werk et al.l I20l4 IStocke et al.l l2013all . Thus, 
understanding the properties of the GGM is critical to 
understanding how galaxies and the IGM co-evolve. 

The low densities in the GGM imply that it is 
very difficult to study in emission. Our ability to 
study the GGM in the present-day universe has been 
dramatically improved with the installation of the 
Gosmic Origins Spectrograph (COS) on the Hubble 
Space Telescope (HST). This has enabled us to use 
the rich suite of absorption-lines due to resonance 
transitions of many important atoms and ions that are 
present in the vacuum ultraviolet. Indeed, over the 
last few years, stu dies such as the COS-Halos, COS- 
Dwar fs , and others (iTumlinson et al.ll2ni laHTriT)D et ahl 
2011t iP rochasl a et al.l 1201 It iB orthak ur et al.l 2013 


Tumlinson et al. 20131: IStocke et ahT 2013a : Werk et alJ 

2014 Bordoloi et al.l I20l4 iLiang fc ChenI I20l4 
Johnson et al.l I2015H have established the connec¬ 
tion between the GGM and many of the key properties 
of galaxies such as star-formation rate (SFR), stellar 
mass, color etc. 

However, so far we have yet to explore the connection 
between the interstellar medium (ISM) and the GGM. 
This is important because in the simplest picture we 
would expect that inco ming gas from th e GGM would 
be largely ionized (e.g. IWerk et aI1l2014ll and upon ac¬ 
cretion would then have to pass through an HI phase and 
then a molecular phase before being converted into stars. 
Thus, the first step in assessing the role of accretion from 
the GGM in the evolution of galaxies would be to look 
at the relationship between the properties of the GGM 
and of the atomic gas in the galaxy disk. 

Of particular interest in this context is a better un¬ 
derstanding of the processes that produce a clear bi¬ 
modality between the population of lower mass galax¬ 
ies undergoing significant star-formation and higher mass 
quiescent galaxies. First discovered at low-redshift (e.g. 


Kauffmann et al.]l2003l iBlanton et al.l[2003l:lBaldrv et ai.l 


2004) . this split into a blue star-forming main sequence 
and a red and dead population is now kno wn to be in 
place out to redshifts of two a nd beyond fe.g. lElbaz et al.1 
I 2 OIII : IWhitaker et al.l l2013h . The galaxy bimodality 
suggests that galaxies undergo fundamental changes 
(through accretion and merging) in their gas content and 
their ability to form new stars after either the galaxy it¬ 
self or its dark matte r halo exceeds a critical mass limit 
('e.g. lLillv et al.l[2013f) . This quenching of star formation 
may be related to a change in the nature of the GGM, 
from a pathway for rapid {vin Vmr) accretion of rel¬ 
atively cold gas (T << Tyiriai) to slow and inefficient 
accretion of hot gas (T ^ Tyir i ai) mediated by radiati ve 
cooling fe.g. iKeres et al.l[2005l . 1200^ iDekel et al.l 1200^ . 

To help probe the processes responsible for the dras¬ 
tic drop in star formation in massive galaxies and the 
implied change in the cold gas supply, we have un- 
dertaken the GALE X Arecibo SDSS Survey (GASS; 
iGatinella et al.1 1201011 . The GASS project combines 
21 cm H I spectroscopic data obtained with the Arecibo 


telescope, optical images and spectroscopy from the 
Sloan Digital Sky Survey (SDSS), and ultraviolet(UV) 
imaging with the Galaxy Evolution Explorer (GALEX) 
to measure the atomic gas and stellar content (star for¬ 
mation rate and history), gas phase metallicity, stel¬ 
lar morphology, and disk rotation speeds of about one 
thousand galaxies over the stellar mass range from 
M* ^ Mq (roughly centered on the transition 

mass between star-forming and quiescent populations). 
We have also obtained molecular ga s data from IRAM 
(GOLD GASS: ISaintonge et al.ll2oilll and long-slit opti¬ 
cal spectroscopy (jMoran et al.ll2012| l for a portion of the 
GASS sample. 

Our sample of galaxies in this paper was drawn from 
the GASS parent sample. This sample enables us to 
probe the connection between the ISM and GGM in L* 
galaxies in a statistically significant manner for the first 
time. In addition, our sample is also complementary 
to the COS-Halos sample as both the samples probe 
galaxies of similar stellar masses. While GOS-Halos 
probed the inner GGM, our sample extends the GGM 
coverage out to the virial radii. Detailed descriptions 
of our sample selection criteria, the GOS observations, 
and data reduction are presented in Section 2. The re¬ 
sults are presented in Section 3. Finally, we summa¬ 
rize our findings and discuss their implications in Sec¬ 
tion 4. The cosmological parameters used in this study 
are Hq = 70 km s~^ Mpc”^, Dm = 0.3, and Ha = 0.7. 

2. OBSERVATIONS 
2.1. Sample: COS-GASS 

Our sample of galaxies was derived from the GASS 
Survey. We have observed all the Quasi Stellar Object 
(QSO) that probe GASS galaxies within a projected sep¬ 
aration from the galaxy (impact parameter) of < 250 kpc. 
The limiting flux of the background QSO was selected 
to be GALEX Far-UV AB magnitude, FUVmag < 19.0. 
This yielded an observed sample of 47 QSO/galaxy pairs. 
The observations of two of those sightlines were found 
to be unsuitable for this analysis. The first is a Broad 
Absorption-Line (BAL) QSO in which it is difficult to 
perform robust measurements of the absorption features 
associated with the target galaxy. The second sightline 
yields a damped Lya system which subsequent obser¬ 
vations show not to be directly associated with the pri¬ 
mary GASS target. We will discuss this unique system in 
Borthakur et al. (in preparation). Hence for the remain¬ 
der of the paper we concentrate on our final GOS-GASS 
sample of 45 QSO sightlines probing the GGM of low-z 
galaxies. 

FigurelUpresents the distribution of our QSO sightlines 
as a function of impact parameter and orientation with 
respect to the target galaxy. The sightlines are marked 
with filled green circles. Those probing face-on galaxies 
being marked in lighter green. The average optical and 
H I disk sizes of the galaxies are shown blue and red 
respectively. The optical sizes for the galaxies were esti¬ 
mated from the SDSS photometric parameter, Rgg. The 
physical size of the Arecibo beam (full width at half max¬ 
imum, FWHM, of 3.3') at the average redshift of 0.037 
for our sample is shown as the yellow circle. The Arecibo 
beam covers the optical extent of every galaxy, but does 
not extend to the virial radius. 
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Fig. 1.— Overview of the COS-GASS sample. The figure shows a representative galaxy from the COS-GASS sample at the center of 
the plot and QSO sightlines shown as filled circles. The sample consists of 45 QSO-galaxy pairs probing the CGM in galaxies from the 
GASS project. The stellar masses of the sample galaxies range between t Mq. The median optical (Rgo) and Hi radii (Rhi) of 

the galaxies in the sample are shown as blue and red ellipses respectively at the center of the plot. The physical size corresponding to the 
FWHM of the Arecibo beam at 21 cm at the average redshift for our sample (z = 0.037) is shown in yellow. The QSO sightlines cover a 
range of impact parameters from 63-231 kpc. The position of the QSO sightlines with respect to the foreground galaxy are shown as filled 
green circles. The orientations of sightlines with respect to the disks are estimated using position angle of the galaxy major axis from SDSS 
photometry. We consider galaxies with de Vaucouleurs a/b parameter greater than 0.7 as having a face-on orientation. These galaxies were 
assigned an orientation parameter of zero degrees and are identified in this figure with a lighter shade of green. 


The principal properties of the galaxies are sum¬ 
marized in Table I. Except where noted below, 
these parameters are taken from the GASS project 
(see ICatinella et al.l l2010L I2012L 120131 for details). 
GASS also incorporat es data products from the SDSS 
(jEisenstein et al.l 120111 ). For example, position and red- 
shift information are based on SDSS photometry and 
spectroscopic measurements respectively. The errors in 
the redshift measurements are of the order of 10“^, which 
corresponds to ~ 30 km s“^ . The values are presented 
in Table 1 and have been rounded off to the forth decimal 
place. The stellar masses of the galaxies in our sample 
range from loiO-i-H-i ]y[Q_ estimated dark matter 
halo masses using the stellar mass mea surements from 
GASS and the methodology described in iBehroozi et al.l 
(j20ir)t) (see Eq.21 and Table 2). The halo mas ses were es- 
timated for z=0 (see for a detailed discussion IShu n iMi 
and they range from 11.6-13.3 Mq. Based on these halo 
masses we estimated the virial radii. Our virial radius 
estimates are within 0.1 dex o f those based on the lu - 
minosity based prescription in iProchaska et al.l (|2011l) . 
We report the halo masses and their corresponding virial 
radii in Table 1. The SFR and by extension the specihc 
SFR (sSFR = SFR/M*) for the GASS sample were de¬ 
rived using a combination of GALEX FUV and NUV 
and SDSS u, g, r, i, z photome try and SDSS spectral-lin e 
indices, (see section 2.2 of ISchiminovich et ahl l2010l ). 
Any value for sSFR < 10“^^ yr~^ can be regarded as 
a upper limit. We quote an uncertainity of 0.2 dex in 
sSFR for the galaxies with sSFR > 10“^^ yr“^. These 
are based on results shown and discussed in detail by 


ISchiminovich et ahl (|2010t) . By including UV fluxes (from 
GALEX) in addition to Ha (from SDSS), the derived 
star formation rates are sensitive to timescales of the or¬ 
der of 10® yrs or larger. These rates define the recent star 
formation activity that might have impacted the GGM. 
For example, star-formation driven outflows traveling at 
200 km s“^ would reach our closest sightline at an impact 
parameter of 63 kpc in about 3x10® yrs and sightlines 
at impact parameter of 200 kpc in 10® yrs. 

Information on the QSO sightlines is presented in Ta¬ 
ble 2. The azimthual orientations of the QSO sightlines 
with respect to target galaxies were calculated from the 
position angle of the optical major axis of the galaxy 
from SDSS r-band photometry. We consider galaxies 
with the de Vaucouleurs AB parameter greater than 0.7 
to be face-on and assign an orientation parameter of 0° 
in the analysis presented in the remainder of the paper. 

2.2. HI masses from single-dish measurements 

Our 21 cm H I spectroscopic data were primarily de¬ 
rived from the GASS observations obtained with the 
single-dish Arecibo telescope. The GASS sample also 
included a few galaxies fr om the Arecibo Le gacy Fast 
ALFA (ALFALFA) Surv ey (iH avnes et al.ll201lD and Gor- 
nell H I archive (GHA, ISoringob et al.l 120051) that met 
the GASS sample selection criteria. GASS HI param¬ 
eters were measured by smoothing the H I profiles fol¬ 
lowed by a baseline subtraction. The edges and peaks 
of the H I profiles were interactively identihed. Fur- 
th er details on th e anal ysis procedures were described 
bv ISoringob et ahl (120051) and Gatinella, Haynes & Gio- 
vanelli 2007 (AJ 134, 334). The velocities at zero flux 
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were marked interactively as part of the width and flux 
measurement process, at or near the position where the 
fits to the edges of the profile cross the baseline. The 
velocity width were corrected for cosmological redshift 
an d instrumental broade ning using procedure described 
in ICatinella et al.l (j20intl . These measurements did not 
include corrections for effects due to inclination or tur¬ 
bulent motion. The measurements of H I masses and 
widths from the GASS database are included in Table 1. 

For five galaxies in our sample, Arecibo H I measure¬ 
ments were not available (although they were members of 
the GASS parent sample). For these we performed simi¬ 
lar observations with the Robert G. Byrd Green Bank 
Telescope (GBT)^. The GBT data were obtained un¬ 
der program GBT-14A-377. We used the dual polariza¬ 
tion L-band system with two intermediate frequency (IF) 
modes and nine- level sampling. The IFs were set to yield 
a channel width of 1.56 kHz (0.33 km s“^ ) using 8192 
channels over a total bandwidth of 12.5 MHz. Observa¬ 
tions were made in the standard position-switching ON- 
OFF scheme with 300 s at each of the positions. Data 
were recorded at 10 s intervals to minimize the effect of 
RFL The OFF position was chosen -|-20' offset in Right 
Ascension from each of the sight lines. During each of 
the observing sessions, we used one of the three standard 
flux calibrators 3G 48 (16.5 Jy), 3G 147 (22.5 Jy), and 
3C 286 (15.0 Jy) - for pointing and estimating antenna 
gain. Local pointing corrections (LPCs) were performed 
using the observing procedure AutoPeak and the correc¬ 
tions were then automatically applied to the data. This 
should result in a pointing accuracy of 3". It is worth 
noting that the pointing accuracy does not signify the 
positional accuracy of the source of the H I flux. The 
H I emitting region may be present anywhere within the 
beam. The flux calibrations were performed by applying 
antenna gain to data for each session separately. Our cal¬ 
ibration error is expected to be no more than 10% . The 
data were analyzed using the NRAO software GBTIDL. 

The GBT observational setup was very similar to that 
of the GASS program and the data products were mea¬ 
sured at the same spectral resolution as the GASS sam¬ 
ple. The H I masses, the velocity widths and the origin of 
the H I data are presented in Table 1. The widths repre¬ 
sent the full velocity extent of the features i.e. the width 
of the 21 cm transition at zero flux. The line-widths were 
estimated visually and the uncertainties associated with 
these values are typically ~20 km s“^ . The H I 21 cm 
velocity ranges cover the systemic velocity of the galaxy 
in all cases. We define the systemic velocity as the op¬ 
tical velocity that traces the stars. The stars represent 
the bulk of the measurable baryons in our galaxies. BY 
doing so, we also alleviates the problem of identifying 
the systemic velocity from H I 21 cm profiles that may 
be asymmetric. 

2.3. Cosmic Origins Spectrograph Observations 

We carried out ultraviolet spectroscopic observa¬ 
tions of the background QSOs with the COS aboard 
the HST. The data were obtained with the grating 
G130M of COS, providing a resolution R= 20,000 — 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 


24,000 (FWHM ^ 14 km s ^ ). The data covered 
observed-frame wavelength of 1140-1470 A. This in¬ 
cluded transitions such as Lyct (A1215.670 A), C II 
(A1334 A), Si II (AI190,1193, &1260 A), Si III (AI206 A), 
Si IV (AA1393,1402 A), and O I (A1302 A). 

The data were reduced using the standard COS 
pipeline. All absorption features in the spectra were 
identified and were exhaustively matched to identify the 
species and transition. Absorption features associated 
with the Milky Way, the target galaxies, intervening sys¬ 
tems, and the background QSOs were identified. We 
use a velocity window of ± 600 km s“^ with respect to 
the systemic velocity (based on SDSS redshifts) to as¬ 
sociate absorption features with our target galaxies. In 
other words, transitions detected within ± 600 km s“^ 
of the systemic velocity of the galaxies are assumed to 
be associated with the target galaxies. The large veloc¬ 
ity window allows us to look for kinematic signatures of 
possible inflowing or outflowing material. Furthermore, 
the velocity window is same as used in the other large 
CGM studies such as COS-Halos (iTurnlinson et al]l2013fl 
and COS-Dwarfs (|Bordoloi et al.1120141 1. thus making our 
analysis consistent with existing CGM studies. 

We do miss a few transitions associated with the target 
galaxies due to blending with other absorption features 
such as from the Milky Way or intervening absorbers be¬ 
tween the QSO and us. In this paper we focus on the 
Lya transition probing the hydrogen in the CGM. A pa¬ 
per describing the metal line properties for the galaxies is 
in preparation. The Lya equivalent width measurements 
are provided in column 9 of Table 2. We fitted Voigt pro¬ 
files to the Lya absorption features to estimate the veloc¬ 
ity centroid of each of the components contributing to the 
full profile. Most of the Lya absorption features (25/36) 
could be described with a single component. There are 
10/36 Lya absorption profiles that required two compo¬ 
nents and one that required three. The velocity centroids 
of all the components (listed in order of the strength of 
the component) and the entire widths of the Lya absorp¬ 
tion profiles within the rest-frame of the galaxy are also 
provided in Table 2. 

3. RESULTS 

We detected Lya absorption in the CGM of 36 out of 
39 galaxies where measurements could be made. In the 
remaining 6 sightlines, the data were contaminated by 
Milky Way or other intervening absorption systems. For 
the three non-detections, we report a limiting equivalent 
width that corresponds to 3 times the noise in the spectra 
in the vicinity of the expected transition. Most of the 
Lya features are saturated (zero residual intensity in the 
line core) with an equivalent width of more than 0.3 A. 

3.1. Tracing the Circumgalactic Medium 

We begin by assessing the connection of the detected 
Lya absorbers to the CGM. Are we seeing gas that is 
physically connected to the galaxy (a true CGM) or sim¬ 
ply material in the inter-galactic medium located close 
to the galaxy in projection? 

First, we can compare the distribution of the Lya ab¬ 
sorbers with respect to the virial radius (Rvir) of the dark 
matter halo. Our sightlines probe gas at impact param¬ 
eters comparable to the Rvir (so, we sample mostly the 









Connecting the CGM to the ISM in galaxies 


5 



Log ( p [kpc]) 




Fig. 2. — Top: Variation of Lyo equivalent widths as a function of impact parameter i.e. p. The COS-GASS data are shown as circles 
and the COS-Halos data are shown as diamonds. The galaxies are divided into two bins based on their specific star-formation rates (see 
text for details). This is indicated by the color of the symbols. Blue and cyan indicate blue star-forming galaxies from the COS-GASS and 
COS-Halos samples respectively. Red and yellow indicate the red passive galaxies. As shown with the black line, the radial distribution 
of the Lyo equivalent width can be fit as exponential with an e-folding scale length of 136 kpc. The grey shaded region represents the 
uncertainty associated with the fit.Bottom: Variation of Lya equivalent widths as a function of normalized impact parameter i.e. p/Rvir. 
The panel on the left presents data using a logarithmic radial scale where as the one on the right is a linear radial scale. The best-fit 
exponential distribution is illustrated with the black solid line. The uncertainty associated with the fit is shown in grey. The best-fit 
parameters are labeled at the bottom of the right panel and imply an e-folding scale length of 1.0 Rvir- 


outer region of the halo). In order to cover the full range 
in impact parameter more completely and with better 
statistic s, we combine our sam ple with the COS-Halos 
sample (|Tumlinson et al.]l2013[) . The COS-Halos sample 
covers a similar range in stellar masses and virial radii 
as our COS-GASS sample. For consistency, we recalcu¬ 
lated the halo masses and virial radii for the COS-Halos 
sa mple from their stella r masses using the prescription 
bv iBehroozi et al.l (|201Clf) . Similar to our limiting values, 
we assign the non-detections in the COS-Halos sample a 
value equivalent to 3 times the root mean squared (rms) 
noise in the spectra. The only remaining difference that 
we do not correct for is the difference in redshift between 
the two samples of about ^0.1. The combined sample 
together covers a large range of impact parameters from 
10-250 kpc. This corresponds to a range in normalized 
impact parameter (p/Rvu) of 0.2-1.4. Thus, the com¬ 
bined sample fully samples the CGM of the galaxies as 
a function of normalized impact parameter. 

We re-investigate the distribution of neutral hydrogen 
as traced by Lya equivalent width as a function of im¬ 
pact parameter and the normalized impact parameter for 
our combined sample. Similar st udies have been done 
in the past with different samples (|Lanzetta et al.l[i995l : 


Chen et al.l 

19981 iTrinn et al.l 119981: IChen et ahl 120011: 

Bowen et al 

120021: iProchaska et ahl 120111: IStocke et ahl 

2013bl: ITumlinson et ahl 120131; iLiane & ChenI 12014 and 


references therein), and the Lya equivalent width was 
found to decrease with increasing impact parameter. Fig¬ 
ure [5] shows the equivalent width of the Lya absorp¬ 
tion line as a function of impact parameter and normal¬ 
ized impact parameter (p/Rvir) for this combined sam¬ 
ple. The Lya equivalent widths are plotted in logarith¬ 
mic units to show the range of values detected. The 
COS-GASS sample is plotted as filled circles and the 
COS-Halos sample is plotted as diamonds. Blue and 
cyan represents galaxies that have specific star formation 
rates (sSFR) of higher than 10“^^ yr“^. These galaxies 
are also referred to as “blue” galaxies in the rest of the 
paper. In red and yellow, are galaxies from the com¬ 
bined sample that have sSFRs lower than that value. 
These galaxies are referred to as “red” g alaxies hereafter. 
The co lor definition is adopted from iTumlinson et ahl 
(l2011bl) . However, the color assignments for the COS- 
GASS sampl e is also consistent with t he more rigorous 
definition bv iSchiminovich et al.l (|2007l see their Fig 7). 
The COS-GASS sample on a plot similar to Fig 7 of 
ISchiminovich et all (|2007ll is shown in Figure [3l The pa- 
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Fig. 3. — Galaxy color definition for the COS-GASS sample. 
The COS-GASS sample on the param eter space of stellar mass 
(M*) and sSFR. Based on the study bv ISchiniinovich et al.l II20071 . 
see their Fig 7), we divide the parameter space into two-zones cor¬ 
responding to red and blue galaxies. The galaxy grouping into the 
two colors f ollowing this defin i tion is consistent with the color def¬ 
inition from ITumlinson et al.l 11201 Ibl i. which defines galaxies with 
sSFR higher than 10“^^ as blue and less than this value as 

red. 

rameter space is divided half-way into red and blue galax¬ 
ies. However, it is worth noting that galaxies closest to 
the dividing line are quite diverse and might correspond 
to a transitional phase (also known as “green valley”). 
For example, some of these galaxies show the Balmer se¬ 
ries transition in emission as well as absorption. They 
also have a wide range of Lya strength in their CGM. 
Unfortunately, our sample does not have a statistically 
significant number of such transitional phase galaxies so 
we can not probe the complexities of their CGM in detail. 

The upper panel of Figure [2] shows that the equivalent 
widths of the Lya absorbers start to drop off beyond an 
impact parameter of about 150 kpc. This same behav¬ 
ior is seen in the plot vs. normalized impact parameter 
(lower left panel). An alternative way to visualize the 
radial distribution of the Lya absorbers is shown in the 
lower right panel, where the normalized impact parame¬ 
ter is plotted in linear units. A linear fit to the data in 
this plot yielded the following best-fit parameters 

LogWpya = (-0.432±0.092)-^-(0.026±0.039) (1) 

rivir 

This implies that the strength of Lya data can be de¬ 
scribed as an exponential function of normalized impact 
parameter with an e-folding length of Rvir as follows: 

WLya = 0.941 e (2) 

This same fit is shown in the lower left panel. A similar 
ht using the (un-normalized) impact parameters is shown 
in the top panel and yields a scale-length of 136 kpc. 

We can also assess the velocity offsets between the 
Lya absorbers and the systemic velocity of the central 
galaxy. Is this gas likely to be bound to the galaxy? We 
have therefore measured the displacement of the centroid 
of the Lya absorption line from the galaxy’s systemic 
velocity. The systemic velocity is defined as the velocity 
corresponding to the SDSS-based redshift measurement 


i.e. Vsys = c z, where c is the speed of light in vac¬ 
uum. The uncertainty in systemic velocity is of the order 
^30 km s“^ . 

Among the sightlines where Lya was detected, 34/36 
(94%) of the absorption features showed the centroid of 
the strongest component of the Lya feature to be within 
250 km s“^ of the galaxy’s systemic redshift. This sug¬ 
gest that the absorbers are most likely bound within 
the dark matter halo of the galaxy. To investigate this 
further, we estimate the escape velocity, for a parcel of 
gas located at the (1) distance of closest impact of the 
sightline (vf^^.) and (2) virial radius (v^"''). The val¬ 
ues are presented in Table 2. These calculations were 
made assuming a s pherically symmetric NFW profile 
(jNavarro et al.l[l996l i for the dark matter halo with a con¬ 
centration parameter of 15. Only 2/36 of Lya absorbers 
have the velocity centroid of the strongest component 
larger than either of the escape velocities. This strongly 
suggests that the bulk of the Lya absorbers are tracing 
the bound GGM. Similarly, if we consider all the compo¬ 
nents of the Lya profiles, then the number of possibly es¬ 
caping components increases to 5/48. In summary, 34/36 
(94%) of the strongest Lya components or 43/48 (90%) 
of all Lya components have velocities smaller than es¬ 
cape velocity. Therefore, statistically the Lya absorbers 
in our sample are likely to be probing the bound GGM. 

We also find that the Lya absorption-lines show veloc¬ 
ities consistent with that of the velocity range observed 
in the 21cm line profiles of the H I gas within the galax¬ 
ies. An example of the 21cm H I profiles of one of our 
galaxies (J1354-I-2433) is presented in the left panel of 
Figure IH The Arecibo spectrum of the galaxy shows 
H I spread over a w ±200 km s“^ (marked as the shaded 
region) in the rest frame of the galaxy. The rest-frame 
velocity was determined by subtracting the systemic ve¬ 
locity {vsys = c z) from the observed H I velocities. The 
velocity centroid of the Lya absorption feature (vuya) 
with respect to the systemic velocity of the host galaxy is 
marked as the red Hlled (strongest component) and open 
(weaker) components. The maximum velocity extent of 
the 21cm H I line profile in the same velocity end as the 
Lya centroid is defined here as vhi. The need to mea¬ 
sure vhi on either the redshifted or blueshifted side of the 
spectra based on the location of the Lya centroid arises 
due to the fact that 21 cm H I profiles may be asymmet¬ 
ric with respect to the systemic velocity. For instance, if 
the Lya centroid is blueshifted in the rest-frame of the 
galaxy, then we measure vhi as the maximum velocity 
of the 21cm H I profile at zero flux blueward of the rest- 
frame. Hence, the estimates of the velocity offsets of the 
ISM and the GGM at the rest-frame are not impacted 
by the asymmetries in the 21cm H I profile. 

The right panel of Figure |4] shows the rest-frame ve¬ 
locity distribution of the Lya absorbers as a function of 
velocity distribution of 21cm H I prohles that traces the 
ISM of the host galaxies. In this figure, the centroids 
of the Lya absorption features are marked with circles 
and the full range of velocity for the absorption features 
are shown as vertical bars. The centroid of the strongest 
component is shown as filled circle and the weaker com¬ 
ponents, if present, are shown as open circles. The color 
of the symbol represents the color of the galaxy as sum¬ 
marized in Table 1. The velocities are normalized to the 
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Fig. 4. — Left: 21 cm Hi spectra of one of our target galaxies (Mhi = 1-8 X 10® Mq, p =155 kpc) from the GASS sample. This 
illustrates a typical Hi profile for our sample. The abscissa and ordinate show the velocity in the rest-frame of the galaxy (based on SDSS 
emission-line estimation) and the fiux in the 21 cm Hi hyperfine transition respectively. The part of the spectrum showing the emission 
(observed, uncorrected for instrumental broadening) is indicated with the grey shaded region (). The spectrum is smoothed to 20 km/s, 
and hence its resolution. The normalized Lya profile for the galaxy is shown in green. The velocity centriod of the strongest component 
of the Lyo; profile is shown as the filled red circle and that of the weaker component is shown as red open circle. The entire width of the 
Lya profile is shown as the red line passing through the circles representing the centroids. The velocity offset of the centroid of the primary 
Lya component from the systemic velocity (at 0 km s“^ in the plot) is termed as v^yo,. We compare this with the velocity extent (half) of 
the Hi profile, vhi-, again with respect to the systemic velocity of the galaxy, towards same velocity direction as the VLya. For example, if 
'^Lya > 0, then we define vhi as the velocity extent of the 21 cm Hi profile towards the positive velocity wing of the profile. For the case 
shown above, Vfjj is defined on the negative velocity wing of the Hi profile and is marked on the plot as the blue double headed arrow. 
Choosing the right velocity edge of the 21 cm Hi profile to compare with the offset of the Lya centroid from the systemic velocity reduces 
the effect of asymmetries on the Hi profile in the analysis. For symmetric 21 cm Hi profiles, v^j is simply the velocity half width of the 
profile. Right: The offset of the velocity centroid of Lya absorbers VLy^ as a function of the velocity (half) widths, vhi, of the 21 cm Hi 
profiles. The velocities are measured with respect to the systemic velocity of the host galaxies based on the redshift measurements from 
optical emission lines. The Hi velocity widths are corrected for instrumental broadening. The centroid of the strongest component is shown 
as a filled circle and those of weaker components are shown as open circles. The full widths of the Lya features are shown as the vertical 
bars passing the centroids. The color of each circle indicates the color of the galaxy. The shaded grey region indicates the parameter space 
where the Lya centroids lie within the Hi extent of the galaxies. The uncertainties in the velocity measurements are shown in the right 
hand corner. 
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Fig. 5. — Left: Lya equivalent width as a function of orientation of the QSO sightline with respect to the optical major axis of the galaxy. 
The open circles represent galaxies that are face-on and their orientation parameter is set to zero. The color of the circles represents the 
color of the galaxy that is summarized in Table 1. Right: The excess in Lya equivalent width ([Log W - Log WjLyo:)? which is independent 
of the impact parameter of the sightlines, as a function of orientation of the QSO sightlines with respect to the target galaxy. The dotted 
line corresponds to the expected equivalent width based on fit to Figure 2. We find no correlations between orientation and Lya strengths. 
The face-on galaxies are excluded while measuring the strength of the correlation. 
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rest-frame using optical redshifts from the SDSS. 

The grey shaded area represents the region of the 
space where |i5VLya| < |5 Vhi|- 90% (26/29) of the cen¬ 
troids of the strongest Lya absorbers lie within the ve¬ 
locity range of the 21 cm H I emission-line associated 
with the target galaxy. Two notable exceptions are the 
blue galaxy J1541-I-2817 and the red galaxy J1515-I-0701. 
The H I spectral range in J1541-I-2817 extends from 
—73 to 107 km s“^ and the Lya absorber was detected 
at -399 km s“^ . Similarly, the Lya velocity centroid of 
J1515-I-0701 is -362 km s“^ whereas the H I associated 
with the galaxy spans a range of —207 to 303 km s“^ 
with respect to the systemic velocity of the galaxy. The 
origin of the offsets in these two systems is not clear 
(they do not seem to be in unusual environments rel¬ 
ative to the other systems). For most galaxies, that is 
20/29 (>66%), the entire Lya prohle lies within the 21cm 
H 1 profile. Therefore,we conclude that the vast majority 
of the absorbing systems are very likely to be gravitation¬ 
ally bound material inside the dark matter halo. 

Next we have looked for any dependence in the proper¬ 
ties of absorbers on the orientation of the sightline with 
respect to the major axis of the galaxy. This would po¬ 
tentially allow us to distinguish between Lya arising in 
an extended disk, in bipolar outflows, or in material sim¬ 
ply tracing the dark matter halo. Figure [S] presents the 
strength of Lya equivalent width as a function of ori- 
entation of the sightl i ne with respect to th e galaxy (e.g 
iKacprzak et al.l[201Cll : iBordoloi et al]l2011h . No correla¬ 
tion has been observed. Since we showed previously that 
Lya equivalent widths correlate with the normalized im¬ 
pact parameters, we corrected for any impact parame¬ 
ter dependence by dividing the observed Lya equivalent 
width by the average Lya equivalent width at that im¬ 
pact parameter (described by the best-fit line of Equa¬ 
tion [T]). We refer to this term as the impact parameter 
corrected excess in equivalent width, or simply as the ex¬ 
cess equivalent width. Here again, we see that the excess 
in the Lya equivalent width does not correlate with ori¬ 
entation of the sightline with respect to the major axis of 
the target galaxy for the full sample. Part of the reason 
for the non-existence of any correlation could be because 
we are probing mostly the outer CGM. 

Therefore, we conclude that the Lya absorbers we 
are studying trace a roughly spherical distribution of 
gravitationally-bound gas with a characteristic size com¬ 
parable to the virial radius. This is what we would expect 
for the CGM. 

3.2. Correlation between properties of the CGM and the 

ISM 

We now explore the connection between the gas in the 
CGM as traced by the Lya absorbers and the cool neu¬ 
tral ISM as traced by the HI 21cm line. Figure [5] shows 
the results for the COS-GASS sample. The abscissa in¬ 
dicates their 21 cm H I mass on the bottom and their 
estimated H I radii on the top. The radius of the disk, 
R-hi was indirectly estimated using the extremely tight 
relationship between H I mass and size as observed by 
ISwaters et al.l (|2002h . We find a positive correlation be¬ 
tween the equivalent width of Lya absorbers in the CGM 

^ Rhi is defined as the annuli of H I mass density of 1 Mq pc“^ 
or column density 1.3 X 10^® cm“^. 


and M(HI). A Kendall r (|Brown et al.lll97'^ test implies 
that the null hypothesis of no correlation can be rejected 
at the 98.5% confidence level^. This also implies a posi¬ 
tive correlation between the equivalent width of Lya and 
the size of the H I disk. We hnd even stronger correla¬ 
tions between the suitably normalized quantities^. This 
is shown in the right panel where the excess Lya equiv¬ 
alent width is plotted as a function of H I mass fraction 
(M(HI)/M*) of the host galaxy. The resulting correla¬ 
tion is inconsistent with the null hypothesis at the 99.8% 
confidence level^. At the same time, we do not find any 
correlation between the Lya equivalent width and the 
stellar mass of the galaxies. 

To visualize our data and to help in its interpreta¬ 
tion, we also performed a stacking experiment to study 
the average properties of the CGM as a function of 
21 cm H I mass within their host galaxies. We di¬ 
vided our sample into four sub-samples based on their 
H I masses. Spectra of each sub-sample were stacked 
within ±1500 km s“^ of the rest-frame Lya transition 
associated with the target galaxy. Each spectrum was 
centered at the centroid of the Lya absorption feature 
before the stacking was performed. Eor the three non¬ 
detections, we used the systemic velocity of the target 
galaxies as zero velocity. The results are presented in 
Eigure[7l Each of these stacks comprise approximately 
8-10 spectra. 

We detected Lya absorption-lines in each stack. How¬ 
ever, their strength shows considerable variations. The 
Lya strengths increased monotonically with increasing 
H I mass within the galaxy. Both the peak depth of the 
absorption features as well as equivalent widths of the 
absorption features increase with increasing H I mass. 
Since most of the individual Lya absorption-lines are sat¬ 
urated, the behavior of the stacks is indicative of the 
increase in the covering fraction of Lya in the GGM 
of galaxies as a function of increasing H I mass. A 
similar trend is also seen as a function of H I fraction 
(M(HI)/M*). Another way to illustrate the same con¬ 
clusion is to note the fractional increase in saturated 
Lya lines as a function of H I mass. Eor example, for 
galaxies with M(HI) < 10®'^ Mq the fraction of Lya ab¬ 
sorbers with WLya > O. 3 A is about 30%. It increases 
to ~55% and 67% for galaxies with M(HI) between 
109.i-9.45 M{m) between Mq respec¬ 

tively. It reaches 100% in the sub-sample with M(HI) 
> 109-9 Mq. 

We also investigated the combined influence of impact 
parameter and H I mass on the equivalent width of the 
Lya profile. To do so, we divided the sample into 2 
impact parameter bins - one lower and the other higher 
than 170 kpc. We then sub-divided each sub-sample into 
galaxies with H I mass less or greater than 9.45 dex. The 
left panel of Figure [8] shows the stacked spectra for each 

® We are using M* as a proxy for halo mass. The normalization 
by the two kinds of masses yield very similar results as Mhaio is 
derived from M*. 

The test was performed on our ce nsored data using the as - 
tronomy survival analysis code ASURV JFeigelson fc Nelsonll 198511 . 
ASURV is capable of handling single and doubly censored data. 
Accuracy of these probabilities can be affected by larger number 
of censored values and other conditions. Since less than a quarter 
of our sample has censored values, we do not expect substantial 
inaccuracies. However, cautio n is appropriate as is for results from 
Kendall’s test on any sample (jWang Welldf20Q(]li . 
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Fig. 6. — Equivalent width of Lyo as a function of 21 cm Hi mass and the excess in Lyo equivalent width ([Log W - Log WjLyct), which 
is independent of impact parameter of the sightline, as a functi on of the Hi mass f raction Mhi/M* of the galaxies. The size of the Hi 
disk is indirectly estimated from the tight Hi size-mass relation ISwaters et aLir2002ll . We find statistically significant correlations with a 
statistical significance of 98.5% for the left panel and 99.8% for the right. The top x-axis of the left panel converts the Hi mass to a radius 
of the Hi disk (see text for details). The average equivalent widths derived from the stacked spectra (presented in the next figure) are 
plotted as green asterisk in the left panel. The averages derived from the data points in the right panel are shown as green asterisk. 
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Fig. 7. — Composite Lya absorption-lines stacked as a function of Hi mass and gas fraction (M(HI)/M*). Individual Lyo absorbers were 
centered at their centroid of the feature i.e. Vcentroid = 0 km s~^ , and then added to produce the stacks. Each of the stacked spectra 
are composed of approximately 8-10 individual spectra in that labeled Hi mass range. The equivalent width and the peak depth (fpeak) 
of the absorption features are labeled in the right corner of each plot in red. The average values of the excess in Lyo equivalent width for 
the sightlines used in creating the stacks are printed in green. These set of stacks are not corrected for impact parameter variations and 
consequently show a weaker trend just like the left panel of Figure 5. Similar plots with impact parameter bins are presented in the next 
figure. 
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Fig. 8. — Top: Stacks of Lya absorption similar to those pre¬ 
sented in Figure [71 but with bins that divide the sample into sub¬ 
sets of impact parameter - one less than and the other greater than 
170 kpc. The values for the bins were chosen such that the sample 
is divided more or less equally into 4 bins. Systematic variations 
the strength of the Lyo absorption as a function of Hi mass within 
each bin can be seen. The equivalent width of the features are 
labeled in the right corner of each plot in red. Bottom: Stacks 
similar to those on the top, but binned by Hi mass fraction and 
normalized impact parameter (p/Rvir)- The sample was divided 
along the median value of normalized impact parameter of 0.70. 
The average values of the excess in Lyo equivalent width for the 
sightlines that go into the stacks are printed in green. As expected, 
these are independent of the normalized impact parameter bins. 

of the bins. This verifies the influence of H I mass as 
well as impact parameter separately on the strength of 
Lya profiles. A similar set of stacks were produced for 
the hybrid parameters - M(HI)/M* and p/Rvir - these 
further confirm the correlation between M(HI)/M* and 
the strength of the Lya absorption lines. 

We also explore the effect of the size of the H I disk 
on the properties of the Lya absorbers. Figure [9] plots 
the equivalent widths of the Lya absorbers as a func¬ 
tion of the impact parameter normalized by the HI ra¬ 
dius (i.e. p/Rhi). Sightlines probing the CGM within 
~7 Rhi all show strong saturated Lya absorption fea¬ 
tures (Wuya > O. 4 A). The other striking feature of the 
plot is the steep drop in the strength of the absorption¬ 


lines at larger normalized radii. 

In conclusion, there is a significant correlation between 
the strength of the Lya absorption-lines tracing the CGM 
and the mass and by extension the extent of the distribu¬ 
tion of 21 cm HI component of the galaxy disk. 

3.3. Relationships of the CGM to the HI and Star 
Formation 

While the properties of the CGM traced by Lya cor¬ 
relate with the H I properties of the host galaxy, the 
H I properties a lso correlate with the star formation rates 
of the galaxies (|Schiminovich et ^120101 : ICatinella et al.l 
l20ir)tl . In this subsection we examine the correlations be¬ 
tween these parameters in an attempt to isolate the pri¬ 
mary physical correlations. This can be addressed most 
easily by determining if there is a correlation between the 
Lya equivalent width and the star-formation rate, and 
then comparing this to the strength of the correlation be¬ 
tween Lya equivalent width and H I mass. This will help 
us probe the physics behind the correlations and isolate 
the driving mechanisms that connect gas at distances of 
hundreds of kpc to the central galaxy. 

We therefore repeat the same exercise for the SFR as 
was done for the H I masses in Figure [H The corre¬ 
sponding plots showing the distribution of Lya equiv¬ 
alent width as a function of SFR and the excess in 
Lya equivalent width as a function of specific SFR (sSFR 
= SFR/M*) are presented in Figure [TUI The Kendall’s 
rank correlation test reveals that the probability of the 
null hypothesis of no correlation between Lya equiva¬ 
lent width and SFR can only be excluded at the 89.9% 
confidence level, while the corresponding probability for 
the correlation between excess Lya equivalent width and 
sSFR is 97.5%. Gomparing these results to the results in 
Figure ini implies that the CGM properties correlate bet¬ 
ter with the properties of the atomic gas than those of the 
current star formation. 

This result may not be surprising. Warm ionized gas 
accreted from the GGM has to pass through the atomic 
phase (traced by H I 21cm) before condensing into molec¬ 
ular phase, which can then facilitate star formation. 
The various processes that connect gas inflow to star- 
formation add intrinsic scatter and therefore, weaken the 
correlation. For example, the efficiency of the conversion 
of gas from the atomic to molecular phase, and of molec¬ 
ular gas into stars, may show substantial variation among 
and within galaxies. The dispersion in the correlation be¬ 
tween Lya equivalent width in the GGM and the SFRs 
within the galaxies may indicate lags in the multi-step 
process of accretion of warm gas via the GGM, conden¬ 
sation to atomic gas, conversion to molecular gas, and 
finally star-formation. 

Nevertheless, we do hnd a signihcant difference in the 
Lya equivalent width distribution between red and blue 
galaxies. The blue galaxies show a uniform Lya strength 
prohle as a function of impact parameter whereas the 
red galaxies show a much larger dispersion at all impact 
parameters including small p/Rvir (Figure 2). We do 
find red g alaxies with strong Lya profiles (similar to that 
found bv lThom et al.l[2012D . but almost half of the red 
galaxies either show weak Lya or even no Lya at our de¬ 
tection limit. To further probe the relationship between 
the GGM and star-formation, a paper detailing the anal¬ 
ysis of the properties of neutral hydrogen and metal lines 
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Fig. 9. — Lyo: equivalent width as a function of the impact pa¬ 
rameter normalized by the radius of the Hi disk (p/Rhi)- The ab¬ 
sorbers are shown as colored filled circles with their colors indicat¬ 
ing the color of the galaxy as defined by their sSFR and presented 
in Table 1. Sightlines within about 7 Rhi produce strong sat¬ 
urated Lya lines. In green we plot the covering fraction of strong 
Lya absorbers with equivalent width larger than 0.32 A (-0.5 dex). 
This corresponds to a column density of 6 x 10^^ cm“^ or larger. 
The right ordinate is labeled to indicate the scale. 

from the combined COS-Halos and COS-GASS surveys 
is in preparation. 

4. SUMMARY & IMPLICATIONS 

The CGM of galaxies acts as the pathway for baryons 
and energy to get in and out of the galaxies. The prop¬ 
erties of the gas in the CGM not only hold clues to the 
evolutionary history of the galaxy but also will influence 
the evolution of the galaxy in the future. After gas from 
the CGM is accreted onto the galaxy, it is expected to 
first pass through a phase as atomic hydrogen located 
primarily in the outer disk. 

Therefore, in order to explore the connection between 
the GGM and the atomic gas reservoir in galactic disks, 
we have undertaken an observational program with the 
COS aboard the HST. This program probed the CGM of 
45 galaxies out to an impact parameter of Ri250 kpc at 
the rest-frame of the galaxies. The sample covers galax¬ 
ies in the stellar mass ranges of i Mq, and as 

such, was designed to probe galaxies spanning the stellar 
mass range where galaxy the galaxy population transi¬ 
tions from lower mass star-forming galaxies that are rich 
in cold interstellar gas to more massive quiescent galaxies 
that contain relatively little cold interstellar gas. 

Our sample was chosen from the GASS project, which 
provided Arecibo H I spectra, GALEX photometry, 
SDSS photometry and spectroscopy and all other related 
data products. For 5 of the galaxies in the sample that 
did not have GASS H I spectra, we obtained H I spectral 
measurements with the GET. 

The COS-GASS sample allowed us to study for the first 
time the connection between the CGM and the ISM of 
galaxies. Based on the analysis presented in this paper, 
we found the following: 

1. We detected Lya absorption in 92% (36/39) of the 
target galaxies for which we have useful data at 
the expected position of the Lya transition. Of 
these detections, 94% (34/36) show Lya absorption 


within ±250 km s ^ of the systemic velocity of their 
host galaxy. 

2. Gombining our sample with the complementary 
GOS Halos sample, we found that the radial dis¬ 
tribution of Lya equivalent widths can be fit as an 
exponential function with an e-folding scale-length 
of Rvir • We defined an excess equivalent width for 
each Lya absorption line with respect to the best- 
fit radial distribution. 

3. The equivalent widths of the Lya absorbers do not 
correlate with orientation of the sightline with re¬ 
spect to the optical major axis of the target galaxy. 

4. These three results imply that Lya absorbers 
are tracing a roughly spherical distribution of 
gravitationally-bound gas with a characteristic size 
similar to the virial radius of the dark matter halo. 
This can indeed be described as a Circum-Galactic 
Medium. 

5. We find positive correlations (at a confidence of 
98.5%) between the equivalent width of Lya and 
Mhi. An even stronger correlation (at a confi¬ 
dence level of 99.8%) is present between the ex¬ 
cess (impact-parameter-corrected) Lya equivalent 
width ([Log W — Log Wjuya) and the H I mass 
fraction (Mhi/M*). 

6 . Similar correlations were also seen in Lya spec¬ 
tra stacked as a function of H I mass and mass 
fraction. The underlying reason for the observed 
increase in Lya absorption strength in the stacks 
with H I mass and mass fraction is due to the 
increase in the covering faction of optically-thick 
neutral gas. For example, the covering fraction of 
Lya absorption is almost 100% in the sub-sample 
with H I mass, M(HI) > 10® ® Mq. This correlation 
was observed for sightlines probing both smaller 
as well as larger impact parameters and is inde¬ 
pendent of the correlation between the equivalent 
width of Lya and impact parameter. 

7. We have used the HI masses to estimate the radii 
of the HI disks. We then found that the equiva¬ 
lent width of the Lya line decreases as the ratio of 
p/R-hi increases. Sightlines probing the gas within 
~ 7 Rhi all show strong Lya absorption features 
(equivalent widths > 0.4A). 

8 . The strength of the Lya absorption is also corre¬ 
lated with the SFR in their host galaxies. How¬ 
ever, we hnd the correlations between the equiva¬ 
lent width of Lya and H I mass and between excess 
Lya equivalent width and HI mass fraction to be 
stronger than the corresponding correlations with 
SFR and SFR/M*. 

These results are consistent with a picture in which 
the reservoir of cold gas in galaxies is fed by accretion of 
gas through/from the CGM. In particular, a process that 
removes or diminishes this CGM reservoir would lead to 
a subsequent drop in the cold gas content of the galaxy 
and hence to a diminished star formation rate. 
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Fig. 10. — Left: Variation in the Lya equivalent width as a function of SFR. The absorbers are shown as colored filled circles with their 
colors indicating the color of the galaxy as defined by their sSFR and presented in Table 1. Right: The excess Lyo equivalent width as 
a function of specific SFR. The correlations have a statistical significance of only 89.9% (left panel) and 97.5% (right panel). These are 
weaker than the corresponding ones in Figure 5. 


The results above are suggestive of a gas accretion his¬ 
tory that is more or less gradual and continuous, at least 
among the Hl-rich galaxies. For a highly episodic pro¬ 
cess we would expect to find only a weak or no correlation 
between the CGM and the cold gas content of galaxies. 
A picture of continuous gas accretion is also consistent 
with the result that material producing strong (optically- 
thick) Lya absorption-lines has a near-unit covering fac¬ 
tor in the CGM of the Hl-rich star-forming galaxies in 
our sample. Interestingly, the Hl-poor quiescent galaxies 
show a wider range in the strength of their Lya absorp¬ 
tion features. Many are weak or even non-detections, 
but we do detect relatively strong absorption feature in 
others. This might indicate a more sporadic accretion 
history for them. 

It is noteworthy that the properties of the CGM traced 
by Lya absorption are more closely correlated with the 
properties of HI than with the properties of the star for¬ 
mation. This may be explained by the fact that the 
H I is mostly in the outer disk and should be more di¬ 
rectly related to the effects of accretion from the CGM. 
In addition, this stronger connection suggests that the 
material probed with Lya absorption is primarily gas 
flowing from the CGM to the galaxy, rather than out¬ 
flowing gas driven by feedback provided by the ongoing 
star-formation. This is further supported by the unifor¬ 
mity in the distribution of the Lya absorbers with respect 
to the galaxy major/minor axis, and the consistency be¬ 
tween the velocity spread of the Lya absorbers and the 
21 cm H I in the galaxy disk. 

Having established a connection between the cool gas 
in the CGM and the HI in the galactic disk, it will be 
important to further probe the interface region between 
the outer disk and the inner CGM to look for clues as to 
how and where the accretion/condensation is happening. 
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TABLE 1 

Description of Galaxy Properties. 


Galaxy 

GASS ID 

RA 

Dec 

z 

Va 

sys 

(km s 

M* 

(Log Mq) 

MLlo 

(Log Mq) 

Rvir 

(kpc) 

Mhi 

(Log Mq) 

AVIii 
(km s ^) 

f*-HI 

(kpc) 

Source 

(Survey) 

sSFR** 
(Log yr-i) 

Golor^ 

J0159+1346 

3936 

29.941 

13.781 

0.0441 

13232 

10.1 

11.6 

163 

9.4 

-126 - 147 

16.2 

GASS 

-9.5 

Elue 

J0808+0512 

19852 

122.068 

5.216 

0.0308 

9227 

10.8 

12.5 

324 

<8.9 

— 

8.4 

GET 

-12.0 

Red 

J0852+0309 

8096 

133.229 

3.152 

0.0345 

10351 

10.3 

11.8 

179 

9.7 

-195 - 143 

22.6 

GASS 

-10.1 

Elue 

J0908+3234 

22391 

137.232 

32.576 

0.0490 

14677 

10.5 

12.0 

212 

<9.3 

— 

14.0 

GET 

<-12.3 

Red 

J0914+0836 

20042 

138.684 

8.601 

0.0468 

14039 

10.0 

11.6 

156 

9.0 

-77 - 108 

10.3 

GASS 

-9.6 

Elue 

J0930+2853 

32907 

142.538 

28.898 

0.0349 

10462 

10.5 

11.9 

208 

9.4 

-190 - 185 

15.0 

GASS 

-10.7 

Elue 

J0931+2632 

53269 

142.817 

26.550 

0.0458 

13718 

11.0 

13.0 

454 

<9.2 

— 

12.0 

GASS 

<-12.6 

Red 

J0936+3204 

33214 

144.101 

32.079 

0.0269 

8073 

10.3 

11.8 

187 

8.7 

-31 - 242 

6.9 

GASS 

-11.7 

Red 

J0937+1658 

55745 

144.292 

16.977 

0.0278 

8328 

10.9 

12.8 

399 

8.9 

-151 - 157 

9.1 

GASS 

-10.3 

Elue 

J0951+3537 

22822 

147.937 

35.622 

0.0270 

8091 

10.6 

12.1 

229 

9.9 

-145 - 240 

31.2 

GASS 

-10.4 

Elue 

J0958+3204 

33737 

149.714 

32.073 

0.0270 

8103 

10.7 

12.3 

269 

10.0 

-412 - -146 

31.6 

GASS 

<-12.7 

Red 

J1002+3238 

33777 

150.711 

32.645 

0.0477 

14300 

10.1 

11.6 

162 

<8.9 

— 

9.1 

GASS 

-11.9 

Red 

J1013+0501 

8634 

153.352 

5.025 

0.0464 

13916 

10.1 

11.6 

165 

9.5 

-187 - 200 

17.2 

GASS 

-10.8 

Elue 

J1032+2112 

55541 

158.196 

21.216 

0.0429 

12852 

10.6 

12.2 

245 

9.9 

-218 - 200 

28.3 

GASS 

-10.1 

Elue 

J1051+1245 

23419 

162.827 

12.757 

0.0400 

11982 

10.4 

11.9 

194 

10.0 

-214 - 108 

33.6 

ALFALFA 

-10.0 

Elue 

J1059+0517 

9109 

164.811 

5.292 

0.0353 

10591 

11.1 

13.3 

563 

<9.0 

— 

9.7 

GET 

-11.9 

Red 

JllOO+1210 

23457 

165.048 

12.171 

0.0354 

10606 

10.1 

11.6 

164 

<8.7 

— 

6.6 

GASS 

-10.7 

Elue 

JllOO+1043 

23477 

165.200 

10.728 

0.0360 

10798 

11.1 

13.3 

601 

9.9 

— 

29.4 

ALFALFA 

-11.0 

Elue 

J1115+0241 

5701 

168.789 

2.699 

0.0442 

13262 

10.7 

12.4 

284 

9.8 

-237 - 227 

27.6 

GASS 

-10.9 

Elue 

J1120+0410 

12452 

170.026 

4.177 

0.0492 

14746 

10.8 

12.6 

332 

<9.1 

— 

10.6 

GASS 

<-12.1 

Red 

J1122+0314 

5872 

170.642 

3.244 

0.0446 

13358 

10.5 

12.0 

222 

<9.2 

— 

12.6 

GET 

-12.0 

Red 

J1127+2657 

48604 

171.943 

26.960 

0.0334 

10001 

10.6 

12.1 

239 

9.8 

-187 - 206 

25.6 

GASS 

-11.0 

Elue 

J1131+1553 

29898 

172.954 

15.897 

0.0364 

10924 

10.2 

11.7 

169 

<9.0 

— 

10.1 

GET 

-12.0 

Red 

J1132+1329 

29871 

173.052 

13.492 

0.0342 

10267 

10.2 

11.7 

168 

9.8 

-256 - 251 

25.0 

ALFALFA 

-9.7 

Elue 

J1142+3013 

48994 

175.575 

30.230 

0.0322 

9665 

10.7 

12.4 

287 

10.3 

-196 - 227 

49.4 

GASS 

-10.4 

Elue 

J1155+2921 

49433 

178.903 

29.351 

0.0458 

13718 

10.5 

11.9 

205 

9.3 

-148 - 162 

13.6 

GASS 

-10.3 

Elue 

J1241+2847 

50550 

190.367 

28.791 

0.0350 

10498 

10.3 

11.8 

180 

9.4 

-125 - 177 

16.4 

GASS 

-10.0 

Elue 

J1251+0551 

13074 

192.894 

5.864 

0.0486 

14572 

10.9 

12.6 

353 

10.2 

-257 - 291 

41.5 

ALFALFA 

-10.4 

Elue 

J1305+0359 

13159 

196.356 

3.992 

0.0437 

13094 

10.4 

11.8 

191 

9.4 

-127 - 262 

15.6 

GASS 

-10.8 

Elue 

J1315+1525 

26936 

198.855 

15.423 

0.0266 

7959 

10.7 

12.4 

294 

<8.9 

— 

9.1 

GASS 

<-12.3 

Red 

J1317+2629 

51025 

199.440 

26.486 

0.0450 

13490 

10.3 

11.7 

176 

9.8 

-133 - 101 

27.6 

ALFALFA 

-10.4 

Elue 

J1325+2714 

51161 

201.345 

27.249 

0.0345 

10336 

10.1 

11.7 

166 

9.5 

-123 - 183 

18.1 

GASS 

-9.8 

Elue 

J1348+2453 

38018 

207.142 

24.891 

0.0297 

8903 

10.1 

11.6 

160 

9.5 

-173 - 173 

17.0 

GASS 

-10.5 

Elue 

J1354+2433 

44856 

208.546 

24.556 

0.0286 

8583 

10.1 

11.6 

158 

9.3 

-160 - 187 

13.3 

GASS 

-11.8 

Red 

J1404+3357 

31172 

211.122 

33.953 

0.0264 

7920 

10.3 

11.8 

181 

9.5 

-193 - 232 

17.2 

GHA 

<-12.3 

Red 

J1406+0154 

7121 

211.678 

1.915 

0.0472 

14153 

10.2 

11.7 

175 

8.9 

-188 - 144 

8.7 

GASS 

-11.8 

Red 

J1427+2629 

45940 

216.954 

26.484 

0.0325 

9734 

10.4 

11.9 

200 

<8.7 

— 

6.7 

GASS 

-12.0 

Red 

J1430+0323 

9615 

217.508 

3.398 

0.0333 

9983 

10.2 

11.7 

166 

9.4 

-133 - 154 

16.2 

GASS 

-11.1 

Red 

J1431+2440 

38198 

217.894 

24.682 

0.0378 

11338 

10.7 

12.2 

256 

9.7 

-211 - 247 

23.2 

GASS 

<-12.7 

Red 

J1454+3050 

42191 

223.516 

30.846 

0.0320 

9584 

10.1 

11.6 

164 

8.8 

-173 - 97 

7.8 

GASS 

-9.8 

Elue 

J1502+0649 

41743 

225.517 

6.823 

0.0462 

13859 

10.5 

11.9 

204 

9.6 

-166 - 186 

20.5 

GASS 

-10.2 

Elue 

J1509+0704 

41869 

227.340 

7.078 

0.0414 

12414 

10.1 

11.7 

166 

9.3 

-133 - 176 

14.0 

GASS 

-9.6 

Elue 

J1515+0701 

42025 

228.781 

7.021 

0.0367 

11005 

10.9 

12.7 

368 

9.5 

-206 - 205 

17.5 

GASS 

-11.9 

Red 

J1541+2813 

28365 

235.344 

28.230 

0.0321 

9620 

10.4 

11.8 

189 

10.1 

-63 - 103 

36.2 

GHA 

-9.6 

Elue 

J1544+2740 

28317 

236.034 

27.673 

0.0316 

9482 

10.1 

11.6 

160 

9.1 

-19 - 191 

10.9 

GASS 

<-12.1 

Red 


^ Systemic velocity, Vsys, is the velocity corresponding to the redshift of the galaxy i.e. Vsys — c z, where c is the speed of light in vacuum. 

Using prescription from Behroozi et al. 2010 (Eq. 21). 

^ The full width of the Hi 21cm line in the rest-frame of the galaxy. The systemic velocity was subtacted from the observed Hi velocity to obtaine the rest-frame velocity. 

Using emperical relationship relating Mhi and Rhi as prescribed bv ISwaters et ^ ||2002D . 

® Values are from f^himinovich et al.l ll2010D : |Catinella et all 1120101) . Values for sSFR < io~^^ yr ^ are regarded as upper limits (discussed in ISchiminovich et ahlf^OlOl) . 
^ Galaxies with sSFR > 10“^^ yr~^ are defined as blue galaxies. Galaxies with sSFR below this value are identified as red galaxies. 
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TABLE 2 

Description of QSO Sightlines. 


QSO 

RAqso 

Decqgo 

ZQSO 

Galaxy 

P 

(kpc) 

p/Rvir 

0a 

WLyo 

(A) 

WLy^/ WLyc 

Liya 

(km s“l) 

Av^ 

Lyo; 

(km S“^) 

yd 

(km s“l) 

^esc,p 

(km s“l) 

V® „ . 
esc,Hvir 

(km s“l) 

J0159+1345 

29.971 

13.765 

0.504 

J0159+1346 

102 

0.6 

64 

1.501±0.023 

2.975 

93,358 

-150 - 

- 380 

147 

178 

150 

J0808+0514 

122.162 

5.244 

0.361 

J0808-I-0512 

215 

0.7 

7 

— 

— 

— 



- 

348 

300 

J0852+0313 

133.247 

3.222 

0.297 

J0852-I-0309 

178 

1.0 

67 

0.113±0.015 

0.322 

50 

0 - 

110 

143 

165 

165 

J0909-I-3236 

137.276 

32.608 

0.809 

J0908+3234 

170 

0.8 

21 

0.090±0.014 

0.213 

101 

50 - 

200 

— 

212 

195 

J0914-I-0837 

138.632 

8.629 

0.649 

J0914-I-0836 

189 

1.2 

69 

0.104±0.020 

0.370 

30 

-50 - 

- 80 

108 

133 

143 

J0930+2848 

142.508 

28.816 

0.487 

J0930+2853 

214 

1.0 

42 

<0.126 

<0.372 

— 



- 

189 

191 

J0931+2628 

142.820 

26.480 

0.778 

J0931+2632 

226 

0.5 

77 

0.114±0.013 

0.199 

200 

150 - 

- 250 

- 

541 

423 

J0936+3207 

144.016 

32.119 

1.150 

J0936+3204 

160 

0.9 

o-t 

<0.113 

<0.280 

— 



- 

183 

172 

J0937+1700 

144.279 

17.006 

0.506 

J0937+1658 

63 

0.2 

64 

0.135±0.021 

0.168 

-204 

-300 - 

- -150 

-151 

669 

370 

J0951+3542 

147.850 

35.714 

0.398 

J0951+3537 

226 

1.0 

3 

0.839±0.014 

2.382 

64 

-90 - 

- 200 

240 

212 

211 

J0959+3203 

149.812 

32.066 

0.564 

J0958+3204 

162 

0.6 

o-t 

0.420±0.016 

0.815 

-238,-144 

-320 - 

- -100 

-412 

298 

249 

J1002+3240 

150.727 

32.678 

0.829 

J1002+3238 

119 

0.7 

40 

<0.063 

<0.140 

— 



- 

166 

149 

J1013+0500 

153.325 

5.009 

0.266 

J1013+0501 

102 

0.6 

o-t 

0.445±0.024 

0.875 

-72 

-180 

- 50 

-187 

179 

151 

J1033+2112 

158.270 

21.204 

0.315 

J1032+2112 

214 

0.9 

45 

0.437±0.033 

1.104 

-1 

-100 

- 95 

-218 

238 

227 

J1051+1247 

162.857 

12.796 

1.281 

J1051+1245 

140 

0.7 

73 

0.781±0.022 

1.706 

-37 

-180 - 

- 140 

-214 

201 

178 

J1059+0519 

164.795 

5.327 

0.754 

J1059+0517 

95 

0.2 

30 

0.270±0.022 

0.340 

6 

-100 - 

- 100 

- 

932 

525 

J1059+1211 

164.984 

12.198 

0.993 

JllOO+1210 

171 

1.0 

61 

0.225±0.014 

0.676 

-61 

-150 

- 0 

- 

148 

151 

JllOO+1046 

165.199 

10.770 

0.422 

JllOO+1043 

108 

0.2 

o-t 

— 

— 

— 



— 

980 

562 

J1115+0237 

168.782 

2.633 

0.567 

J1115+0241 

209 

0.7 

83 

0.195±0.020 

0.430 

2 

-100 - 

- 100 

227 

293 

263 

J1120+0413 

170.021 

4.223 

0.545 

J1120+0410 

162 

0.5 

78 

0.830±0.018 

1.434 

201,370 

50 - 

390 

— 

396 

308 

J1122-I-0318 

170.601 

3.301 

0.475 

J1122-I-0314 

221 

1.0 

of 

0.110±0.018 

0.315 

67 

0 - 

150 

— 

205 

204 

J1127-I-2654 

171.902 

26.914 

0.379 

J1127+2657 

140 

0.6 

26 

0.705±0.021 

1.342 

-28,-178 

-250 - 

- 100 

-187 

267 

220 

J113H-1556 

172.905 

15.946 

0.183 

J1131+1553 

176 

1.0 

of 

— 

— 

— 



— 

153 

155 

J1132-I-1335 

173.044 

13.586 

0.201 

J1132+1329 

230 

1.4 

5 

0.319±0.017 

1.318 

57 

-40 - 

- 130 

251 

138 

155 

J1142-1-3016 

175.551 

30.270 

0.481 

J1142-1-3013 

104 

0.4 

50 

0.886±0.023 

1.351 

3 

-200 - 

- 170 

227 

376 

265 

J1155-I-2922 

178.970 

29.377 

0.520 

J1155-I-2921 

208 

1.0 

1 

0.742±0.023 

2.156 

72,-150 

-200 - 

- 230 

162 

188 

189 

J1241-I-2852 

190.374 

28.870 

0.589 

J1241-I-2847 

198 

1.1 

40 

0.211±0.020 

0.667 

33 

-120 - 

- 170 

177 

160 

166 

J1251+0554 

192.853 

5.906 

1.377 

J1251+0551 

200 

0.6 

57 

0.409±0.021 

0.765 

80 

-20 - 

- 180 

291 

401 

328 

J1305+0357 

196.351 

3.959 

0.545 

J1305+0359 

103 

0.5 

11 

0.821±0.016 

1.494 

67,-43 

-160 - 

- 180 

262 

218 

175 

J1315+1525 

198.938 

15.432 

0.448 

J1315+1525 

155 

0.5 

17 

0.405±0.018 

0.728 

64 

-50 - 

- 170 

- 

341 

272 

J1318+2628 

199.508 

26.475 

1.234 

J1317+2629 

198 

1.1 

86 

0.184±0.032 

0.597 

0 

-120 - 

- 120 

-133 

155 

162 

J1325+2717 

201.266 

27.289 

0.522 

J1325+2714 

199 

1.2 

55 

— 

— 

— 



— 

143 

153 

J1348+2456 

207.093 

24.947 

0.293 

J1348+2453 

153 

1.0 

81 

0.474±0.035 

1.301 

-97 

-230 

- 0 

-173 

150 

147 

J1354+2430 

208.604 

24.502 

1.878 

J1354+2433 

155 

1.0 

78 

0.545±0.034 

1.537 

-97,-5 

-170 

- 50 

-160 

147 

146 

J1404+3353 

211.118 

33.895 

0.549 

J1404+3357 

111 

0.6 

57 

0.749±0.027 

1.469 

-26 

-150 - 

- 150 

-193 

198 

166 

J1406+0157 

211.732 

1.954 

0.427 

J1406+0154 

222 

1.3 

67 

- 

- 

- 



- 

147 

161 

J1427+2632 

216.898 

26.537 

0.364 

J1427+2629 

170 

0.9 

of 

- 

- 

- 



- 

195 

184 

J1429+0321 

217.420 

3.357 

0.253 

J1430+0323 

231 

1.4 

of 

0.807±0.027 

3.414 

-39,109 

-150 - 

- 250 

-133 

135 

153 

J1431+2442 

217.858 

24.706 

0.407 

J1431+2440 

110 

0.4 

18 

0.569±0.015 

0.929 

73,156 

0 - 

220 

247 

317 

236 

J1454+3046 

223.601 

30.783 

0.465 

J1454+3050 

223 

1.4 

37 

0.472±0.035 

1.941 

57 

-50 - 

- 160 

97 

134 

151 

J1502+0645 

225.517 

6.754 

0.288 

J1502+0649 

224 

1.1 

80 

0.438±0.013 

1.384 

12,-55 

-150 - 

- 110 

186 

182 

188 

J1509+0702 

227.368 

7.043 

0.418 

J1509+0704 

130 

0.8 

62 

0.956±0.022 

2.220 

49,-18,-215 

-275 - 

- 130 

176 

167 

153 

J1515+0657 

228.781 

6.952 

0.268 

J1515+0701 

180 

0.5 

14 

0.270±0.023 

0.468 

-367 

-500 - 

- -300 

-206 

439 

342 

J1541+2817 

235.340 

28.285 

0.376 

J1541+2813 

128 

0.7 

of 

0.864±0.011 

1.797 

-363 

-520 - 

- -250 

-63 

200 

174 

J1544+2743 

236.114 

27.723 

0.163 

J1544+2740 

196 

1.2 

55 

0.191±0.022 

0.687 

126 

50 - 

210 

191 

136 

147 


Orientation of the QSO sightlines with respect to the disk of the galaxies. The values are based on SDSS r-band photometric measurements. 

^ Centroid of the multiple components of the Lya absorption feature as estimated via Voigt profile fitting in the order of the strength of the component. 

^ Full width of the absorption feature in the rest-frame of the galaxy. 

Maximum velocity of the Hi profile in the rest-frame of the galaxy. This value was measured on the velocity side where the centroid of the Lyman <y absorber was detected. 
® Line of sight escape velocity at the impact parameter probed by the QSO sightline assuming a NFW profile for the galaxy’s dark matter distribution. 

^ Face-on galaxies. 
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